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Abstract 



We analyze a quantum kinetic equation describing both boson and fermion pair production. We 
explore the solution of the kinetic equation in its Markovian limit. The numerical study shows 
^ ' an enhancement (bosons) or a suppression (fermions) of the pair creation rate according to their 

, quantum statistics. The modification of the time evolution of the distribution function is small 

CO ' but for strong fields more pronounced than for weak fields. 
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The pre-equilibrium evolution of the quark-gluon plasma, believed to be created in an ultrarel- 
ativistic heavy-ion collision, is described by means of a transport equation that can be used to 
Oh! explore dynamical properties in the plasma phase. The formation of the QGP is assumed to pro- 

^ • ceed via the creation of a strong chromoelectric field in the region between the two receding nuclei 

, after the collision. The field subsequently decays by emitting quark-antiquark pairs according to 

i-J^ ' the nonperturbative tunneling process of the Schwinger mechanism. 

The process of parton pair production within the Schwinger mechanism has been addressed 
. . ■ by many authors in recent years with the back reaction scenario also considered. While 

rS I these studies have been very useful in exploring new phenomena such as plasma oscillations, the 

■ phenomenological approaches suffer from the lack of a derivation of a source term within a kinetic 

theory. Indeed such a consistent field theoretical treatment leads to a modified source term pro- 
viding a non-Markovian evolution of the distribution function. This result was first obtained by 
Rau Q for a constant electric field. A generalized treatment allowing for a time dependent field 
was demonstrated in Ref. H. A systematic numerical study of the time structure of such a source 
term was provided in Ref. M for the case of the creation of boson pairs. In Ref. §, using a field 
theoretical treatment we have derived a kinetic equation characterized by the following properties: 
(i) it has non-Markovian character, (ii) the distribution function of the produced particles has a 
non-trivial momentum dependence and (iii) due to the statistical factor [1 ± 2f±{P,t)] the pro- 
duction is modified according to the quantum statistics of the created pairs. The properties of the 
source term itself such as the momentum dependence and the time structure have been studied in 
Ref. |5j for a constant and a time dependent field. However, a numerical analysis of the evolution 
of the distribution function was not provided. In this paper we explore the solution of the kinetic 
equation in the Markovian limit in order to compare fermion and boson pair production. 
Using a simple field-theoretical model of charged fermions in an external, homogeneous, time- 
dependent field characterized by the vector potential A^^ = {0,0,0, A{t)) with A{t) = A3{t) and 
the resulting electric field E{t) = E^{t) = —A{t) = —dA{t)/dt, we obtained the kinetic equation 
for the distribution function f± ~ f±{P,t) 



df± , 9/± e^E{t)sl r , , E{t') 



dt ' 'dP\\{t) 2tuHp,t) uHp,t') 



l±2/±(t') 



cos[x{t' ,t)], (1) 



where the upper sign (lower sign) corresponds to fermion (boson) pair creation. We define the 
total energy uj'^{p,t) = e\ + Puit), the transverse energy £\ = m? -\- p\ and P\\{t) = p\\ — eA{t). 




Furthermore x{t',t) = 2[8(p, t) — Q{p,t')] denotes the difference of the dynamical phases which 
are defined as ^ 

e(p,i)= / dt'iu{p,t') . (2) 

Jto 

The lower limit in Eq. (||), < = to, describes the system in the asymptotic field free state, E{to) = 0. 
The Markovian limit of Eq. (|I|) is defined by the neglect of retardation effects in the source term 
and reads 

^I±^ = [l±2U{P,t)]J{P,t), (3) 
where J{P,t) is the source term independent of the statistics (low density limit) 

Taking into account the initial condition lim f± — 0, we obtain the following solution of the 

t — ^ — oo 

kinetic equation 



/±(i^,t) = T^(l-exp' 



±2 / dt'j(P,t') 



(5) 



This result is exact in the Markovian limit and holds for any time-dependent homogeneous electric 
field. One can show analytically that the distribution functions are positive for all times and 
momenta. In the lowest order of the expansion of Eq. (^) we obtain 

For the numerical estimation we choose the simple case of a constant electric field defined as 
A{t) = tEo. In Fig. 1, we show the source term in Markovian limit (J(P,t)) as function of 
momentum (l.h.s.) and time (r.h.s.). Resulting from the Ansatz of a constant external field, we 
observe oscillations which can be identified as Airy functions In the left panels of Fig. 2, the 
distribution functions (Eq. (||)) for fermions and bosons are plotted for a strong and a weak electric 
field as function of the dimensionless time {te±) for vanishing parallel momenta. Compared to 
/o of Eq. (||), the distribution function for bosons (/+) is enhanced while that for fermions (/_) 
is suppressed. For both fermions and bosons the evolution approaches a constant value at large 
times. For weak fields the difference between /+ and /_ is very small, see lower panel. 
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Figure 2: On the l.h.s. the time evolution of the distribution functions of bosons (/-|_) and 
fermions (/-), Eq. @ compared with /o of Eq. (P) at p\\ = for a strong field (upper 
panels), EQ/e\ = 1.5, and a weak field (lower panels), EQ/e\ = 0.5, is shown. On the 
r.h.s. the corresponding modified production rates are plotted. 



We have analytically and numerically explored the solution of a quantum kinetic equation describ- 
ing particle production in the Markovian limit. We observe that the resulting production rate 
depends on the statistical character of the produced particles. This effect is elucidated in the right 
panels of Fig. 2 where the time evolution of the modified production rate is plotted. We obtain a 
suppression of the production rate for the creation of fermion pairs and an enhancement for boson 
pairs. Although the effect is surprisingly small, it is much more pronounced for strong fields than 
for weak fields. 
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